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Abstract 

We present numerical results on bubble profiles, nucleation rates and time evolution for a weakly 
first-order quark-hadron phase transition in diflterent expansion scenarios. We confirm the stan- 
dard picture of a cosmological first-order phase transition, in which the phase transition is 
entirely dominated by nucleation. We also show that, even for expansion rates much lower than 
those expected in heavy-ion collisions nucleation is very unlikely, indicating that the main phase 
conversion mechanism is spinodal decomposition. 
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1. Introduction 

Although Lattice QCD results seem to indicate a crossover for the quark-hadron tran- 
sition [1], the possibility of a weakly first-order transition is not ruled out yet from the 
experimental point of view and some observables behave differently for a first-order tran- 
sition and for a crossover [2]. Moreover, Lattice QCD does not provide any information 
on the (nonequilibrium) dynamical nature of the transition, whose critical behavior could 
be very different from what one might expect from a crossover in the (equilibrium) phase 
diagram. 

We calculate exact critical bubble profiles and their critical radii and surface tension as 
functions of the temperature [3] using realistic equations of state (EoS). We also compute 
the fraction of the plasma hadronized via nucleation and its temperature as functions of 
time in order to analize quantitatively the importance of nucleation in the dynamics of 
the phase transition. 
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2. A model for the phase transition of a QGP in expansion 



After a hot plasma of quarks and gluons (QGP) is formed, it cools down due to 
expansion and eventually reaches the critical temperature T,.- From this point on, the 
plasma supercools and, as it becomes metastable, nucleation starts to act effectively. 
When one critical bubble is actually formed and the corresponding region of the plasma 
is hadronized, the latent heat released in the process may reheat the plasma. This goes 
on until the whole plasma is converted to the hadronic phase and then the temperature 
suddenly falls. However, if the expansion is fast enough, the released heat cannot keep the 
system from supercooling down to some temperature under which it becomes unstable 
and suffers spinodal decomposition. In this case, usually called a quench, nucleation is a 
secondary mechanism in the phase conversion [4|5j . 

Our criterion for determining the importance of nucleation relies on the fraction /(t) 
of the plasma which is hadronized from the moment T = Tc and before the unstable 
spinodal temperature is reached. 

We assume local thermal equilibrium and (approximate) entropy conservation. For the 
quark phase {Nf = 2 + 1), we use two different equations of state: one corresponding to 
lattice results [6] and the other one to a Bag Model [7]. For the hadronic phase, we take 
a gas of over 250 free massive resonances [S] . We use the 3-d Euclidean action [S] 



= j d\ 



i[V0(r)]2 + y[0(r),T] 



(1) 



where is a phenomenological effective potential given by 

V{(j), T) = a{T)<j)^ - bTf + (2) 

We assume that the quark phase corresponds to the order parameter 4'q{T) = and 
the hadron phase to (j)h{T). The connection with thermodynamics is made through the 
identification pq(T) —pj(T) = V{(j)j,T) for a homogeneous phase with </> — (pj. 

We also suppose a homogeneous and isotropic (Hubble-like) expansion of the plasma, 
with scale factor a{t) = aoexp{HQt), where Hq is a constant. A final assumption is the 
homogeneous reheating of the expanding plasma due to the release of latent heat. All 
these ingredients together lead to the equation [lOj 



m = 



a{t) J ■'^ ' s{T{t)) ^ ' 



(3) 



where s = (1 — f)sq + fsh is the spatial average of the entropy density, Sq(T) (sh(T)) is 
the entropy density of the quark (hadron) phase and As{T) = Sq{T) — Sh{T). 

3. Method for calculating T(t) 

3.1. Bubble features 

Critical bubbles are spherically symmetric extrema of the action ([Ij, so that 

^^l^^W ^'^"^ m = ^h and cj^{oo) = cbq=0. (4) 
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This equation may be solved analytically for temperatures close to Tc in the so-called 
thin- wall approximation [719] . However, this approximation, which is valid for i?c/^ 2> 1 
~ T~^), fails for temperatures much lower than Tc and one has to solve eq. (|4]) 
numerically. This is what we do in order to calculate the main ingredients for the dynamics 
of the phase conversion: the critical bubble radius Rc{T) and surface tension cr(T). 

3.2. Dynamical quantities 

In order to calculate /(t), a necessary ingredient for T{t) according to eq. we use 

m 

f{t)^l-expl^-jdt' (^^Jr[T{t')]^R'{t',t)^, (5) 
where the nucleation rate per unit volume per unit time is [7] 

and the radius of a bubble that was born at time t' with radius Ro[T{t')] and grew until 
time t with velocity v^iT) = -V{(l)h,T)/[aiT)T/2] is [10] 

t 

Rit',t) = ^fo[T{t')] + J dt" v^[T{t")]^^. (7) 

t' 

4. Results and discussion 

We have solved numerically the set of equations ([3]) - ([7]) for the temperature as a 
function of time: T(t). 
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Fig. 1. T{t) for different values of the expansion 
rate Hq. 



Fig. 2. T(i) for different equations of state 
(H-^ = 600fm/c). 



3 



In Fig.m T{t) is shown for three different expansion rates: very fast = lOOfm/c), 

intermediate {Hq^ = 600fm/c) and very slow {Hf^^ — lO^fm/c). Notice that the plasma 
suffers a quench into the spinodal region in the fast expansion, for it does not have enough 
time either to nucleate a non negligible number of bubbles or to reheat. In the very slow 
case, the system presents a very slight supercooling which is quickly washed away by the 
plasma reheating, so that the phase conversion follows quite closely the equilibrium, as 
expected from a very slowly expanding system. Inverse expansion rates of the order of 
600fm/c (for the Hadron+Lattice, or HL, EoS) lead to the most interesting scenarios, 
in which the plasma reaches a considerable supercooling but reheats before reaching the 
spinodal temperature. After this reheating, the temperature once again approaches Tc 
and the phase conversion is completed mainly by the percolation of the bubbles born in 
the supercooling stage. 

In Fig. [2] we show the effect of the latent heat on the time evolution of the temperature. 
The HL EoS leads to a stronger supercooling and a shorter coexistence time in comparison 
with the Hadron+Bag (HB) EoS, which has a larger latent heat. Once the HL EoS leads 
to a smaller latent heat, it takes longer to reheat the plasma, and it spends more time in 
a low temperature (and high nucleation rate) regime. This means that when the system 
is successfully reheated, a large fraction of it is already hadronized and the remaining 
fraction is not able to maintain the expanding system with T <Tc for much time. On the 
other hand, when the latent heat is larger (e.g., with the HB EoS), the plasma reheats 
sooner and can keep its temperature high for a longer time. 

5. Conclusion 

Using realistic equations of state (in comparison with the usual Bag Model EoS), 
we showed that if QCD has a weakly first-order phase transition, then nucleation is 
very unlikely the main phase conversion mechanism in HIC experiments, regardless of 
reheating effects. In this scenario, the adequate description for the transition should be 
a quench followed by spinodal decomposition. In the early Universe, however, we showed 
that nucleation completely dominates the conversion dynamics. Finally, we showed that 
a stronger transition leads to a smaller supercooling and to a larger coexistence time. 
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